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FLOW-STRUCTURE  INTERACTION  NOISE 
AT  LOW  MACH  NUMBERS 

M.  S.  Howe 

Boston  University,  College  of  Engineering 
110  Cummington  Street.  Boston  MA  02215 


SUMMARY 

The  method  of  compact  Green’s  functions  is  described  for  calculating  the  sound  produced 
by  low  Mach  number  turbulence  interacting  with  a surface  and  by  surfaces  in  relative  motion. 
The  method  is  very  powerful  and  can  supply  analytical  predictions  of  complex  fluid-structure 
interaction  noise,  and  can  also  be  used  to  make  acoustic  predictions  by  ‘postprocessing’  data 
obtained  from  a numerical  simulation  of  a hydrodynamic  interaction.  The  theory  is  illustrated 
by  cases  of  sound  production  by  a vortex  interacting  with  a sphere,  by  the  linear  and  nonlinear 
theories  of  parallel  vortex-airfoil  interaction  noise,  the  trailing  edge  noise  produced  by  an  airfoil  of 
finite  thickness  and  of  arbitrary  chord,  and  by  the  compression  wave  generated  when  a projectile 
is  fired  into  a long  duct. 


1.  EQUATION  OF  VORTEX  SOUND 

Lighthill’s  acoustic  analogy  [1]  can  be  recast  to  emphasize  the  importance  of  vorticity  (u»)  and 
entropy  fluctuations  as  sources  of  sound  by  taking  the  total  enthalpy  B to  be  the  fundamental 
acoustic  variable  [2,  3].  Lighthill  identified  the  Reynolds  stress  as  the  principal  acoustic  source  in 
homentropic  flow.  In  that  case 
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where  p is  fluid  density,  p = p(p)  the  pressure,  and  v denotes  velocity,  and  Lighthill’s  equation 

KprOTTlP‘? 

<L1> 

where  c is  the  local  speed  of  sound. 

In  the  irrotational  acoustic  far  field  Crocco’s  form  of  the  momentum  equation  dv/dt  — -VB 
implies  that  B = -d( p/dt,  where  y?(x,  t ) is  the  velocity  potential  that  determines  the  whole  motion 
in  the  irrotational  regions  of  the  fluid.  B is  therefore  constant  in  steady  irrotational  flow,  and  at 
large  distances  from  the  sources  perturbations  in  B represent  outgoing  sound  waves.  If  the  mean 
flow  is  at  rest  in  the  far  field,  the  acoustic  pressure  p = p0B  = -p0d<p/dt. 

At  low  Mach  numbers  M local  mean  values  of  p and  c differ  from  their  uniform  respective  values 
p0  and  c0  by  terms  of  relative  order  M2.  The  vortex  sound  equation  (1.1)  can  then  be  simplified 
(a)  by  taking  c = c0  and  p = p0,  and  (b)  by  neglecting  nonlinear  effects  of  propagation  and  the 
scattering  of  sound  by  the  vorticity.  The  production  of  sound  is  then  governed  by 

(^“v2)B  = div(l,,Av)'  M 


where  in  the  acoustic  field 

P(x,t)  ~p0B(x,t).  (1.3) 
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2.  VORTEX-SURFACE  INTERACTION  NOISE 

Consider  the  solution  of  (1.2)  in  the  presence  of  solid  boundaries  S that  are  either  stationary  or 
vibrating  at  small  amplitude.  Introduce  a Green’s  function  G(x,y,t  - r)  that  satisfies 

/ 1 ft  ft  \ 

[ft0Wft~dyi)G  = S{x- y)S{t-T)'  G = 0forT>*’  (2-D 

and  has  vanishing  normal  derivatives  dG/dxn,  dG/dyn  respectively  for  x and  y on  S.  Then 
Green’s  theorem  [3  - 5]  supplies  the  solution  of  (1.2)  in  the  form 

/dG  t dG 

(w  A v)(y,  r)  • — (x,  y,  t - r)d3ydr  + v ^ u>(y,  t)  A — (x, y,  t - r)  • dS(y)dr 

r Q 

+ js  G(x,  y , t - r)  (y,  r)dSj  (y )dr,  (2.2) 

provided  frictional  effects  on  S are  effectively  the  same  as  for  incompressible  flow,  where  v is  the 
kinematic  viscosity,  and  the  vector  surface  element  dS  is  directed  into  the  fluid.  The  first  integral 
describes  the  production  of  sound  by  convecting  vorticity.  Green’s  function  takes  account  of  the 
influence  of  the  body  on  the  efficiency  of  sound  production  by  this  source.  The  second,  surface 
integral  represents  the  sound  produced  by  frictional  forces  on  the  body;  it  can  be  neglected  in 
most  applications  at  high  Reynolds  numbers.  The  final  term  is  the  contribution  from  normal 
velocity  fluctuations  on  S. 

Therefore,  at  large  Reynolds  number  and  for  a non-vibrating  surface,  the  pressure  in  the 
acoustic  far  field  is  given  by 

p(x,  t)  = -p0  |(uAv)  (y,  r)  ~ (x,  y,  t - r)d3ydr.  (2.3) 

The  integral  includes  contributions  from  the  direct  radiation  by  quadrupole  sources  within  the 
fluid  (with  intensity  oc  p0v3M 5 per  unit  volume  [1])  and  also  from  dipole  sources  on  S.  For  a 
body  with  surface  irregularities  comparable  in  size  to  the  turbulence  length  scale  (when  the 
irregularity  is  acoustically  compact)  the  dipole  intensity  oc  p0u3M3;  for  turbulence  near  the  leading 
or  trailing  edge  of  an  airfoil  of  non-compact  chord  the  intensity  typically  scales  like  p0v3M2. 
Thus  at  low  Mach  numbers,  compact  surface  irregularities  (or  a compact  foreign  body)  and 
edges  of  non-compact  bodies  are  usually  the  dominant  sources  of  flow  generated  sound.  We  shall 
consider  these  cases  separately  and  demonstrate  how  detailed  acoustic  predictions  can  be  made  for 
complex  interactions  in  terms  of  the  appropriate  compact  approximation  for  the  Green’s  function  G. 

3.  COMPACT  BODIES  AND  COMPACT  SURFACE  IRREGULARITIES 

3.1  Compact  bodies  and  cylindrical  bodies  of  compact  cross-section  The  leading  order 
dipole  radiation  is  calculated  by  approximating  G in  (2.3)  by  the  compact  Green’s  function  [2,  3] 

C(x'y'‘-T)  = (31) 

X = x-<p*(x)  \ 

> Kirchhoff  vectors  for  the  body, 

Y = y-V(y)  J 
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where  the  vector  components  Xj(x)  and  Yj( y)  are  the  velocity  potentials  of  incompressible  flow 
past  the  body  having  unit  speed  in  the  j-direction  at  large  distances  from  the  body;  ip*  is  the 
velocity  potential  of  the  incompressible  flow  that  would  be  produced  by  rigid  body  motion  of  S 
at  unit  speed  in  the  j-direction.  For  a cylindrical  body  of  compact  cross-section  parallel  to  the 
^-direction,  we  take 

A's  = £3,  Yi  = J/3  - 

The  approximation  (3.1)  represents  the  solution  of  (2.1)  correct  to  dipole  order  provided  either 
x or  y (or  both)  lie  in  the  acoustic  far  field  of  S,  in  the  limit  in  which  the  acoustic  wavelength 
A > t ~ the  length  scale  of  S (Figure  1).  Special  cases  referenced  in  the  illustrations  given  below 
are  listed  in  Table  1: 


Standard  special  cases 


body 

*1 

*2 

*3 

sphere  of  radius  a 
center  at  origin 

Xl  (1  + w) 

12  ( 1 + W) 

l3(1  + w) 

circular  cylinder  of  radius  a 

~ ( 1 1 “2 

r (1  | a2  'i 

Z3 

coaxial  with  the  ^3-axis 

strip  airfoil 

-a  < X\  < a , x2  = 0,  -oo  < m3  < 00 

II 

Re(— i\/z2  — a2) 
z = x 1 + ix  2 

X3 

Table  1 


■ y 


Figure  1 


Figure  2 - 


3.2  Surface  irregularity  This  case  involves  turbulence  sources  adjacent  to  a surface  that  is 
locally  plane  on  the  scale  of  the  acoustic  wavelength  except  for  the  presence  of  one  or  more 
acoustically  compact  surface  irregularities  (Figure  2).  The  Green’s  function  incorporates  ‘images’ 
of  all  dipole  sources  ‘in’  the  wall,  such  that  the  net  normal  dipole  strength  vanishes  identically.  If 
the  wall  coincides  with  x2  = 0,  with  fluid  in  x2  > 0,  then 


G(x,y,t  -r) 


4tt|X  - Yj 


(- 


IX  — Yh  1 
c0  ; + 47r|X-Y| 


lX-Yj\ 

Co  ) 


1 


(3.2) 


(SYA)  14-4 


where 


Yi  = yi  -<tf(y),  Y2  = y2 , Y3  = y3  - <p*3(y) 
Xi=xi~  v?t(x),  X2  = rr2,  X3  = x3  - (p*3{x)  ^ 


(3.3) 


and  X = (Xi,  — X2,X3).  The  case  of  a two-dimensional  surface  irregularity  (projection  or  cavity) 
uniform  in,  say,  the  ^-direction  is  obtained  by  setting  Y3  = 1/3,  X3  = x3. 


3.3.  Radiation  from  a stationary  compact  body:  vortex  interacting  with  a sphere  To 
evaluate  (2.3)  take  the  coordinate  origin  within  S and  consider  an  observer  at  x in  the  acoustic  far 
field  (|x|  ->  00).  Expand  the  compact  Green’s  function  G to  first  (dipole)  order  in  the  retarded 
time  x • Y/c0|x|: 


G(x,y,t-r) 


47rjX-Y|<5 


r - 


)X-Y|\ 

Co  ) 


4?r|x 


where  the  prime  denotes  differentiation  with  respect  to  t.  The  first  5-function  is  independent  of  y 
and  makes  no  contribution  to  (2.3),  which  becomes 


M-*00  (3'5) 

The  integral  represents  the  retarded  dipole  strength  produced  by  normal  stresses  on  S.  Indeed,  the 
component  Fj  of  the  unsteady  force  exerted  by  the  body  on  the  fluid  is  given  (for  incompressible 
flow)  by  [3] 

Fj ; = -p0  J u}  A v • VY)d3 y - r?  j>  VYj  ■ u>  A dS,  77  = p0u.  (3.6) 

where  the  second,  frictional  component  of  force  can  be  discarded  at  very  large  Reynolds  numbers. 

To  illustrate  the  procedure  consider  the  sound  produced  when  a line  vortex  of  strength  T 
convects  past  a fixed,  rigid  sphere.  Let  the  sphere  have  radius  a with  center  at  the  coordinate 
origin,  and  let  the  vortex  be  initially  far  upstream  of  the  sphere  and  parallel  to  the  13-axis  at  a 
distance  h above  the  plane  of  symmetry  x2  = 0 (Figure  3).  There  is  a mean  irrotational  flow  in 
the  xi -direction  at  speed  U,  which  at  sufficiently  low  Mach  number  ( M = U/c0)  is  given  by 

U = UVXi(x),  (3.7) 


where  A”i(x)  is  given  in  the  first  row  of  Table  1.  The  vortex  is  convected  towards  the  sphere  by 
the  mean  flow  and  deformed  as  it  passes  around  the  sphere;  more  distant  parts  of  the  vortex  (at 
|x3|  a)  are  unaffected  and  remain  parallel  to  the  ^-direction.  Figure  3 illustrates  the  distorted 
vortex  at  several  values  of  T = Ut/a , when  h/a  = 0.2  and  when  the  self-induced  motion  of  the 
vortex  is  neglected. 

When  surface  viscous  effects  and  vortex  shedding  are  ignored,  the  sound  generated  during  the 
interaction  is  given  by  (3.5)  in  which  u is  approximated  by  the  vorticity  of  the  distorted  line  vortex. 
There  is  no  contribution  from  unsteady  ‘drag’  ( j = 1)  because  v = f/VYi  and  10  A VY)  • VYj  = 0. 
Similarly,  there  can  be  no  net  ‘side-force’  on  the  sphere  because  of  the  symmetric  shape  of  the 
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xi 


Figure  3 


Figure  4 


deformed  vortex,  and  therefore  there  will  be  no  contribution  from  j = 3.  The  sound  is  accordingly 
produced  by  a dipole  orientated  in  the  ^-direction;  i.e.  the  interaction  produces  an  unsteady  ‘lift’ 
in  this  direction  which  is  responsible  for  the  sound,  and  which  can  therefore  be  cast  in  the  form 


p(x,t) 


-p0U  cos  0 d 
47tc0|x|  dt 


|x|  -» oo, 


(3.8) 


where  0 = cos 


(x2/|x|)  is  the  angle  between  the  observer  direction  x and  the  X2-axis. 


Typical  plots  of  the  nondimensional  pressure 


P(x,  t) 

p0n7Mcos©/47r|x| 

are  shown  plotted  against  U[t]/a  in  Figure  4 for  two  values  of  h/a ; they  illustrate  how  the  sound 
level  decreases  as  the  initial  standoff  distance  h of  the  vortex  increases. 

In  this  calculation  the  unsteady  interaction  has  been  evaluated  in  a linearized  approximation, 
when  the  unsteady  drag  is  predicted  to  vanish.  Thus,  if  the  sphere  is  ‘stuck’  to  a plane  wall  in  the 
presence  of  a turbulent  flow,  the  present  approximation  would  yield  no  ‘dipole’  sound,  since  only 
the  drag  and  side  force  can  then  make  a finite  dipole  contribution  because  (from  (3.2)) 


(?(x,y,f-T) 


'27t|x| 


C0/  27TCo|x|J  V Co 


(3.9) 


3.4.  Parallel  blade-vortex  interactions  The  same  method  of  solution  in  terms  of  the  compact 
Green’s  function  is  applicable  to  airfoils  of  compact  chord  and  arbitrary  span,  but  it  is  sufficient 
here  to  give  explicit  results  when  the  span  is  compact. 
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Figure  5 


Consider  a planar  airfoil  of  rectangular  or  elliptic  planform,  orientated  as  illustrated  in  Figure  5 
at  zero  angle  of  attack  to  a mean  flow  at  speed  U in  the  x\ -direction.  A spanwise  line  vortex  of 
strength  F is  swept  past  the  airfoil  at  an  initial  standoff  distance  h above  the  airfoil,  as  indicated 
in  the  side  view.  For  an  airfoil  of  compact  chord  the  acoustic  pressure  produced  by  the  interaction 
is  given  by  equation  (3.5)  where,  for  an  airfoil  of  span  L (between  -\L  <x 3 < ~L), 


Yi=y  1,  Y2 


Re(-iv/^  — a(y3)2),  |y3|  < \L 

2/2,  M > 


Yz  = y3,  z = yi+m, 


where  2a  (y3)  is  the  airfoil  chord  at  the  spanwise  location  3/3.  For  the  rectangular  airfoil  0(1/3)  = a = 
constant;  for  the  elliptic  airfoil  we  set 


«(lft) 

a 


N < \L- 


(3.10) 


Vorticity  will  be  shed  into  the  wake  of  the  airfoil  in  accordance  with  the  Kutta  condition.  This 
smooths  out  conditions  at  the  trailing  edge,  so  that  sound  is  generated  primarily  by  the  interaction 
of  the  vortex  with  the  leading  edge.  This  is  accounted  for  by  removing  the  trailing  edge  singularity 
of  the  Green’s  function,  by  modifying  the  X2-component  of  Y as  follows: 

Y2  = Re(]/2a(y3)]/z  + a(y3)'),  |y3|  < ^L.  (3.11) 


In  the  linearized  approximation  v = (17,0,0)  in  (3.5),  and  the  line  vortex  remains  rectilinear; 
the  radiation  is  produced  entirely  by  the  lift  dipole.  If  the  vortex  crosses  the  midchord  of  the 
airfoil  at  time  t = 0,  (3.5)  yields 


p(x,t) 

p0TUM  cos  0(L/|x|)/4jra 


(3.12) 
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where  yz  = yz/L.  For  a rectangular  airfoil  the  right  hand  side  evaluates  to 


— rim 
25 


(*?+! +%y) 


The  acoustic  pressure  signatures  (the  left  hand  side  of  (3.12))  for  the  rectangular  and  elliptic 
airfoils  are  plotted  in  Figure  6 for  a vortex  standoff  distance  h/a  = 0.2.  For  the  rectangular  airfoil, 
the  peak  amplitude  is  larger  and  the  width  of  the  acoustic  pulse  is  narrower  than  for  the  elliptic 
airfoil,,  where  the  interaction  of  the  leading  edge  with  the  vortex  occurs  over  a longer  time  period. 


U[t]la 

Figure  6 


When  the  h — 0 ‘image’  vorticity  in  the  airfoil  prevents  the  vortex  from  impinging  on  the 
leading  edge,  and  causes  the  trajectory  to  be  locally  deflected  above  the  airfoil  (for  T > 0).  For  a 
rectangular  airfoil  this  case  can  be  treated  by  assuming  that  the  section  of  the  vortex  within  the 
span  —\L  < £3  <\L  follows  a path  that  is  approximately  the  same  as  for  locally  two-dimensional 
flow,  and  by  neglecting  small  contributions  from  the  ends  of  the  airfoil.  There  are  now  two  distinct 
dipole  components  of  the  sound:  that  discussed  above  ( P2 , say)  associated  with  the  unsteady  lift, 
and  a second  smaller  component  (p\)  produced  by  the  unsteady  ‘suction’  force’  at  the  leading 
edge.  The  corresponding  nondimensional  pressures 

Pi  fat) P2fat) 

p0YUM  cos<3>(L/|x|)/47ra’  p0TUM  cos©(£/|x|)/47ra 

are  plotted  in  Figure  7 for  a velocity  ratio  T/4i:aU  = 0.2.  The  upper  part  of  the  figure  shows  the 
path  followed  by  those  sections  of  the  vortex  inboard  of  the  airfoil  tips. 
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4.  NONCOMPACT  BODIES 

4.1  Edge  noise  Our  discussion  will  be  framed  in  terms  of  the  canonical,  parallel-sided  airfoil 
shape  illustrated  schematically  in  Figure  8,  which  has  chord  t and  uniform  thickness  h,  and  a 
- rounded  nose.  Take  the  origin  0 at  a convenient  point  near  the  trailing  edge  such  that  the  ‘upper’ 
and  ‘lower’  planar  surfaces  of  the  airfoil  are  at  £2  = ±5/1,  the  a; i -axis  is  in  the  direction  of  the 
mean  flow,  and  £3  is  parallel  to  the  airfoil  span  (out  of  the  plane  of  the  paper  in  Figure  8a). 

When  these  conditions  are  satisfied  the  far  field  acoustic  pressure  fluctuations  p(x.u)e~'ut 
of  frequency  u>  produced  by  the  interaction  of  the  turbulence  with  the  airfoil  are  given  by  the 
following  form  of  (3.5) 

/dG  \ 

— (x,y,o>)- (u/Av)(y ,w)d3y>  (4-1) 


where  G(x,y,u;)  is  the  time  harmonic  Green’s  function,  which  satisfies  G(x,  y.t  — r)  = 


At  low  Mach  numbers  the  trailing  edge  noise  is  dominated  by  contributions  to  (4.1)  from 
the  neighborhood  of  the  edge.  The  wavelength  of  the  sound  ~ S/M,  where  5 is  a length  that 
characterizes  the  eddy  size  near  the  edge,  and  may  in  practice  be  of  the  order  of  the  boundary 
layer  thickness  or  the  airfoil  thickness  h.  In  either  case  the  acoustic  wavelength  greatly  exceeds 
h when  M«l,  which  implies  that  the  turbulence  responsible  for  the  edge  noise  is  always  very 
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much  closer  than  an  acoustic  wavelength  from  the  edge.  Then  the  nondimensional  source  distance 
K0\Jy\  -f  y\  ~ \Jy\Jr  /acoustic  wavelength  <C  1,  where  k0  = w/c0  is  the  acoustic  wavenumber. 
When  the  observation  point  x is  in  the  acoustic  far  field  and  y is  within  the  source  region,  Green’s 
function  can  be  expanded  in  terms  of  this  small  parameter.  When  the  airfoil  chord  £ is  noncompact 
(k0£  1),  the  first  nontrivial  term  in  this  expansion  yields 


G(x,y,w) 


-Y^sina  ^sin  </j*{y)e”s‘>lx  5,3,3 
ny/2id\x  - y3i3j 


k0£  ->  oo,  K0\fy(+yi  < 1.  (4.2) 


where  i3  is  a unit  vector  parallel  to  the  airfoil  span  (in  the  positive  23-direction).  Mr  — + x\ 

denotes  the  observer  distance  from  the  edge  of  the  airfoil  (~  the  23-axis),  then  (r,9)  are  the  polar 
coordinates  of  the  observer  position  relative  to  a plane  23  = constant,  i.e.  (21,22)  = r(cos0,sin0); 
ip  = sin-1(r/|x|)  is  the  angle  between  the  observer  direction  x and  the  edge  (see  Figure  8b).  The 
function  ip*( y)  = <p*{yi,  1/2)  of  the  source  position  y is  the  velocity  potential  of  ideal,  incompressible 
flow  around  the  edge  (in  the  anticlockwise  direction),  and  therefore  depends  on  the  geometrical 
shape  of  the  edge.  At  source  distances  yjy\  + y|  from  the  edge  that  are  large  compared  to  the 
airfoil  thickness  h but  small  compared  to  the  chord  £,  <p*  (y)  must  tend  to  the  following  expression 
for  the  potential  of  flow  around  a rigid  half-plane : 

(p*(y) Vrism(8'/2),  h<£r'<g.£,  where  {y  1,2/2 ) = r'(cos0', sin 9').  (4.3) 


In  applications  where  the  wavelength  is  not  necessarily  small  compared  to  the  chord  l the 
representation  (4.2)  is  modified  in  its  dependence  on  frequency  and  on  the  directivity,  but  the 
edge  potential  function  ip*( y)  is  unchanged.  This  is  illustrated  in  Figure  9 [6]  for  ip  — 90°  and 
k0£  =1,  5,  10  and  50.  The  heavy  curves  are  polar  plots  of  the  linear  pressure  amplitude 
(normalized  to  the  same  maximum  value);  in  each  case  the  light  curve  corresponds  to  the  sin  ) 
directionality  for  k0£  = 00.  At  the  lowest  frequency  shown  in  Figure  9 {k0£  — 1)  the  radiation 
peaks  in  directions  normal  to  the  airfoil,  and  the  directivity  coincides  with  that  for  a dipole  source 
orientated  in  the  22-direction  (as  for  the  compact  chord  airfoil  of  §3.4).  Multiple  lobes  develop 
as  the  frequency  increases,  and  the  directivity  tends  in  an  oscillatory  manner  towards  that  for 
the  semi-infinite  airfoil.  In  all  cases  there  are  radiation  nulls  in  the  directions  0 = 0°  and  180°, 
respectively  downstream  and  upstream  of  the  airfoil. 

The  above  argument  implies  that  the  influence  of  finite  chord  on  the  trailing  edge  noise  frequency 
spectrum  does  not  depend  on  the  edge  geometry.  In  Figure  10  we  compare  the  predicted  edge 
noise  spectra  for  a semi-infinite  airfoil  with  that  for  an  airfoil  of  finite  chord  l for  two  different 
radiation  directions,  and  for  the  following  parameter  values  that  are  typical  of  model  scale  tests 
conducted  in  air  and  in  associated  numerical  simulations  [7  - 13] 

-f-«78,  M = 0.09,  (4.4) 

<5* 


where  <5*  is  the  boundary  layer  displacement  thickness  in  the  vicinity  of  the  trailing  edge. 
Significant  departures  from  the  half-plane  prediction  occur  for  u5*/U  < 0.4;  in  particular  spectral 
levels  are  much  reduced  at  frequencies  below  that  of  the  spectral  peak,  which  is  shifted  to  higher 
frequencies. 
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4.2  Projectile  entering  a duct  The  overall  pressure  rise  across  the  compression  wave  produced 
ahead  of  a bullet  shaped  body  projected  into  a long,  but  finite  duct  of  cross-sectional  area  A = irR? 
fitted  with  a long,  cylindrical  entrance  ‘hood’  of  length  £h  and  cross-section  Ah  = nR\  (Figure  11) 
is  given  approximately  by  [14] 


A p = 


p0U 2 Ap 

(1-M2)  A 


(4.5) 


where  U and  A0  respectively  denote  the  projection  speed  and  the  cross-sectional  area  of  the  body, 
and  M = U/c0. 

When  the  blockage  A0/ A is  small  and  the  projection  Mach  number  M does  not  exceed  about 
0.4,  the  initial  form  of  the  compression  wave  can  be  calculated  from  the  following  modified  form 
of  (1.2): 


(|  & - ' v2) B = (' + £) I (t£(* + WWr2 w*3)) + div(“ A v)-  (4  6) 

The  new  source  term  is  a slender  body  representation  of  the  body,  which  is  assumed  to  travel 
along  the  axis  of  the  duct  (the  negative  Zi-axis),  where  .4t(s)  is  the  body  cross-sectional  area  at 
distance  s from  the  tip  of  the  nose,  which  is  taken  to  cross  the  entrance  plane  of  the  hood  at  time 
t = 0. 

The  characteristic  wavelength  of  the  sound  produced  when  the  body  interacts  with  the  duct 
entrance  is  sufficiently  large  (~  Rh/M) that  only  plane  waves  can  propagate  into  the  duct.  The 
corresponding  compact  Green’s  function  for  an  observer  at  x deep  inside  the  duct  then  has  the 
form 


G(x,y;f-r 

c0  7j 
2 Ah 


t {h  (w  - r - g*±pB)  + RbH  (h - r - ^-gM) } 


x\  -¥  -oo,  (4.7) 

where  [t]  = t + (zi  - i')/c0  is  the  retarded  time  (f  ~ 0.61  Rh  being  the  ‘end  correction’  of  the  hood 
entrance  [15]).  The  velocity  potential  <p*(x)  represents  an  irrotational,  incompressible  flow  out  of 
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the  hood  from  aq  = -oo,  normalized  such  that 

9?*(x)  ~ xi  — £'  when  |a:i|  jR/j,  within  the  hood, 

-Ah 

~ ■■  ,■■■,  when  |x|  » Rh  outside  the  hood.  (4.8) 

47t|x| 


Also, 


Re  — — 1,  Rj  = 


Ah  — A 
Ah  + A 


2Ah 
Ah  + A 


(4.9) 


Re  is  the  reflection  coefficient  for  long  sound  waves  at  the  open  end;  Rj  is  the  reflection  coefficient 
at  the  junction  of  the  hood  and  uniform  duct  for  reflection  back  into  the  hood;  and  7j  is  the 
transmission  coefficient  for  transmission  from  the  hood  into  the  duct.  Formula  (4.7)  is  uniformly 
valid  for  source  positions  y lying  within  the  hood  and  the  neighborhood  of  its  ends. 


h 


(b) 


displaced  air  ► 


exit  flow 


vorticity 


projectile 


2 h 


VWR 


Figure  12 


Figure  11 

To  understand  the  terms  in  the  brace  brackets  of  (4.7),  consider  the  case  where  the  source  point 
y lies  in  the  body  of  the  hood.  For  n — 0 the  first  Heaviside  function  represents  the  front  of  the 
disturbance  radiated  into  the  duct;  the  second  is  the  contribution  from  the  wavefront  that  initially 
propagates  towards  the  hood  entrance,  where  it  is  reflected  with  reflection  coefficient  Re  = — 1 and 
subsequently  transmitted  into  the  duct;  the  amplitude  of  each  of  these  ‘direct’  waves  is  increased 
on  transmission  through  the  junction  into  the  duct  because  the  transmission  coefficient  7j  > 1- 
The  terms  in  n > 1 represent  contributions  arriving  at  later  times,  after  the  transmission  into  the 
duct  of  those  components  of  the  direct  waves  that  have  suffered  n reflections  from  both  ends  of 
the  hood;  the  amplitude  decreases  by  a factor  Rj  < 1 on  each  reflection  back  into  the  hood  from 
the  junction  with  the  duct,  so  that  the  higher  order  modes  rapidly  decrease  in  amplitude. 

The  amplitude  and  profile  of  the  sound  wave  radiated  from  the  far  end  of  the  duct  when 
the  compression  wave  arrives  is  proportional  to  the  compression  wave  pressure  ‘gradient’  dp/dt 
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incident  on  the  exit.  The  linear  theory  prediction  of  dp/dt  is  given  by  the  solution  of  (4.6)  in 
terms  of  the  Green’s  function  (4.7).  In  a long  duct  it  would  also  be  necessary  to  account  for 
nonlinear  steepening  of  the  compression  wave;  the  solution  of  (4.6)  should  therefore  be  regarded 
as  providing  the  initial  conditions  for  such  a nonlinear  calculation. 

The  principal  source  of  the  compression  wave  is  the  monopole  representing  displacement  of  fluid 
by  the  nose  of  the  projectile  and  a dipole  generated  by  the  drag  produced  by  the  pressure  rise  over 
the  nose,  both  of  which  are  accounted  for  by  the  first  term  on  the  right  of  (4.6)  [14].  The  vortex 
source  represents  secondary  effects  of  separation  on  the  body,  on  the  duct  walls,  and  in  the  shear 
layers  of  the  ‘exit’  flows  from  the  hood  mouth  and  through  the  junction  from  the  duct  into  the 
hood,  of  the  fluid  displaced  by  the  entering  projectile.  The  contribution  from  the  monopole  and 
dipole  can  be  written 


dp 

dt 


PoU 3 
(1  - M2) 


M)IL 


00  dAr  d7G 

T(yi  + Ur)^(xyyuOAt 


dy\ 


dyf 


r)dyidr. 


(4.10) 


Figure  12  illustrates  a comparison  of  this  solution  with  an  experiment  in  air  using  a 7m  long 
horizontal,  unflanged,  circular  cylindrical  duct  made  of  hard  vinyl  chloride,  with  inner  radius  R = 
5 cm.  An  axisymmetric  body  is  projected  into  the  duct,  guided  by  a 5 mm  diameter  taut  steel  wire 
extending  along  the  duct  axis.  Pressure  measurements  are  made  using  a wall-mounted  transducer 
1.5  m from  the  hood  entrance.  The  hood  consists  of  a circular  cylindrical,  vinyl  chloride  pipe  of 
inner  radius  Rh  = 6.25  cm;  therefore  Ah/ A = R\/R?  — 1.5625.  A 10  cm  collar  at  one  end  of  the 
hood  facilitates  an  airtight  and  smooth  mating  with  the  duct.  When  in  place  the  hood  has  an 
overall  length  of  60.3  cm,  and  a ‘working  length’  of  4 — 50  cm. 

The  projectile  had  an  ellipsoidal  nose  profile  obtained  by  rotating  the  curve  y = 
hyJ{x/L)(2  -x/L),  0 < x < L about  the  x-axis,  so  that 


At{x) 

Ao 


!(2-f),  0 < x < L, 
1,  x > L, 


where  h — 2.235 cm  L — 11.18cm, 


(4.11) 


and  A0  — rh7.  The  model  had  overall  length  equal  to  124.3  cm.  and  fore-aft  symmetry  with  equal 
nose  and  tail  profiles 

The  solid  curve  in  Figure  12  is  the  prediction  of  equations  (4.7),  (4.10)  for  U = 294 kph 
(M  = 0.24),  plotted  as  a function  of  the  non-dimensional  retarded  time  U[t]/R,  where  the  nose  of 
the  projectile  crosses  the  entrancejdane  of  the  hood  at  t = 0.  The  open  triangles  are  observed 
values  of  the  pressure  gradient  measured  in  the  duct  at  a distance  of  1.5  m from  the  entrance  plane 
of  the  hood.  The  agreement  between  theory  and  experiment  is  excellent;  the  wavy  nature  of  the 
pressure  gradient  profile  is  caused  by  interference  between  waves  multiply  reflected  (temporarily 
trapped)  within  the  hood  before  transmission  into  the  duct.  The  theory  based  on  (4.7)  is  seen  to 
capture  this  interference  perfectly. 
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